Minimizing driver's irritation at a roadblock 
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^"5 Abstract. Urban traffic is a logistic issue which can have many societal implications, 

especially when, due to a too high density of cars, the network of streets of a city 
becomes blocked, and consequently, pedestrians, bicycles, and cars start sharing the 
same traffic conditions potentially leading to high irritations (of people) and therefore 
to chaos. In this paper we focus our attention on a simple scenario: We model the 
driver's irritation induced by the presence of a roadblock. As a natural generalization, 

(j we extend the model for the two one-way crossroads traffic presented in [5] to that 

of a roadblock. Our discrete model defines and minimizes the total waiting time. 

^ The novelty lies in introducing the (total) driver's irritation and its minimization. 

C/3 Finally, we apply our model to a real-world situation: rush hour traffic in Hillegom, 

The Netherlands. We observe that minimizing the total waiting time and minimizing 

C/} the total driver's irritation lead to different traffic light strategies. 

>■> 

PACS. 45.70.Vn Granular models of complex systems; traffic flow - 05. 40. -a Fluc- 
Oh tuation phenomena - 82.20.Mj Nonequilibrium kinetics 

T—H 

^ 1 Introduction pie highway systems) but are mathematically 
{N( NP complete (and therefore untractable) at the 
Traffic flow problems are a major topic of re- level of cities with thousands of traffic lights, 
search. Examples include traffic on highways or In this paper, we propose a somewhat differ- 
within cities (crossings, traffic light controls, road- ent view and consider the driver's perspective, 
blocks and so on). Many practical questions are We develop a model for the driver's irritation^ 
still in search for an optimal answer: Which traf- and propose to minimize this instead of mini- 
fic light settings lead to a minimal waiting time mizing the total waiting time (TWT). We will 
per car? Can one control this in an a priori see that such a strategy may lead to slightly dif- 

^ way? What about the waiting time for all the 



cars together stacked in a lane? At which speed i There are all sorts of causes for irritation in traf- 

do incoming vehicles have to drive to maximize fie, often actions of other drivers like tailgating, lack 

the flow during a traffic jam? How can local of flashing when switching lanes, driving too slow or 

car-car and car-traffic light interactions lead to too fast and so on. It may also be caused by having 

emergent coordination of patterns like "green- to wait for a long time at traffic lights or in jams, 

waves" ? Most of these questions can be addressed Here we start from the assumption that high irrita- 

for not too complex road networks (like sim- tions lead with an increased probability to sudden 

dangerous motions and herewith to accidents. 
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ferent clearance policies than when minimizing 
the TWT. In next papers, we intend to apply 
this modeling framework to optimizing for load 
levels within capacity and will study the depen- 
dence of optimization on modeling the concept 
of "driver's irritation". 

Within this paper we study one specific traf- 
fic situation: a two-way street roadblock where 
one lane is blocked. Our model is a general- 
ization of that for the intersection of two ur- 
ban streets presented in [5]. Note that there 
are standard measures taken by the authorities 
for a roadblock situation. An example is a by- 
pass: traffic is redirected so that there is less or 
no traffic that has to pass the roadblock. Also, 
roadwork can be planned during time intervals 
when there is not much traffic, such as week- 
ends or nights. Sometimes these measures are 
not possible and one has to look for alternatives 
such as lowering the irritation with the help 
of smart traffic light settings. Typically traffic 
lights are installed such that the TWT is mini- 
mized. This means that the total flow of traffic 
is optimized. 

This work is based on pi. More research on 
traffic flow is done by many researchers, for in- 
stance [3] connects the current research trends 
in traffic flow to discrete crowd dynamics. An- 
other approach is queueing theory, see e.g. [2[7]. 
We use a discrete deterministic particle systems- 
like model pi. A more stochastic approach of 
traffic flow problems is found in [IT] and refer- 
ences cited therein. Finally traffic flow can be 
modeled by fluid-dynamics principles using or- 
dinary and partial differential equations; see, for 
instance, Refs. [XJ|9j . 

In Section [2| we present our roadblock model 
based on |5 . We minimize the total waiting time. 
Here we introduce our concept of driver's irri- 
tation and study traffic light settings that mini- 
mize total irritation. In Section [3j we apply our 
methodology to a real traffic situation, using 
empirical data for rush hour scenario in Hil- 
legom, The Netherlands. Finally, we conclude 
with a discussion in Section [4J 

2 Roadblock model 

Consider a two-way street where one of the two 
lanes is blocked. See Figure [I] for a schematic 
picture of the roadblock. Traffic lights are in- 
stalled such that both directions do not have 



green light at the same time. After each time 
a direction has had green light, both lights are 
red for a specific time. This is called a transit 
period and is necessary to clear the road. This 
leads to the following cycle: 

— Phase I: direction 1 green: period T\ 

— Phase II: both directions red: period r 

— Phase III: direction 2 green: period T — T\ 

— Phase IV: both directions red: period r 

After Phase IV, the cycle repeats itself and starts 
with Phase I again. We take the values for T and 
r constant, and the value for T\ variable. We are 
now looking for a value for T\ that minimizes the 
total waiting time. We note that we can see this 
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Fig. 1. Schematic picture of the roadblock 

problem as a queuing system with two queues 
and one server, where at every moment no more 
than one queue is served. 

2.1 Notation. Basic assumptions 

— The arrival rates 0.1,0.2 [cars/second] and 
the passing rates /?i,/?2 [cars/second] are as- 
sumed to be constant; 

— The cycle number is given by n, and repre- 
sents in which cycle we are; 

— The queue lengths of each direction at the 
beginning are 0. At time t the queue lengths 
are: JVi(t),JV 2 (t); 

— The total waiting time TWT is defined as 
the sum of the waiting times of each direc- 
tion. This waiting time is defined as the prod- 
uct of the total number of cars that have to 
wait and the time that they have to wait 
(times of the red phases); 

— We only consider the heavy traffic scenario 
where during one cycle more cars arrive than 
leave in both directions. This assumption is 
made because this scenario is the most inter- 
esting and problematic one for drivers, logis- 
tics companies and governmental and envi- 
ronmental institutes. Note however, that in 
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this scenario the queue of waiting cars would 
grow infinitely. In reality, the heavy traffic 
scenario holds only for limited time (such as 
during rush hour), which ensures the queue 
to be bounded. 

Here we assume the cycle period T is constant. 
This is a fair assumption for the rush hour case. 
Choosing T\ fixes T 2 , the cycle period T playing 
the role of a parameter. 



2.2 Total waiting time 

Based on these assumptions, we can define the 
total waiting times for Directions 1 and 2 as: 

T^CTi) := Ni(nT+Ti)(T—Ti)+^ai(T—Ti) 2 . 

(1) 

T { n w ' 2 \T x ) :=N 2 (nT)(T 1 + r) + N 2 {(n + 1)T - r 
+ ±a 2 (T 1 +2T) 2 . 

(2) 

The assumption of heavy traffic means that the 
queue lengths after every phase increase linearly 
over the cycle number. For example the queue 
length of Direction 1 is: 

N^nT + T,) := n( ai T - ft^) + (a, - frfr. 

(3) 

Let us substitute these queue lengths into (2). 
Now if we take the sum of those two waiting 
times, we get the total waiting time, that is the 
waiting time of all the cars at the roadblock. 
Minimizing the total waiting time with respect 
to T\ gives the value: 

rr, (ai-a2 + /3i+/3 2 )T-4/3 2 r 

Tl = wrwi • (4) 

We note that if we take r = 0, we get the 
same value as the one obtained in [5] , which sug- 
gests that the roadblock problem and the two 
one-way crossroads have the same mathemati- 
cal structure. 



2.3 Driver's irritation 

Instead of minimizing the total waiting time, we 
can take a different look at the traffic situation. 
We like now to think in terms of minimizing the 



driver's irritation. First of all, we have to define 
our concept of irritation. We say that the to- 
tal irritation per cycle of one direction depends 
on the irritation per car. This irritation per car 
is defined at the moment that the lights switch 
from green to red. We call this the crucial mo- 
ment because at that time there is a significant 
change in the irritation: it increases because the 
people have to wait while they first were able 
to pass the roadblock. We note that the crucial 
moment for Direction 1 is after Phase I and for 
Direction 2 after Phase III. In this framework, 
we assume the change in the irritation per car 
depends on several causes: 

(i) Waiting time; the longer you have to wait 
after the crucial moment, the higher the 
irritation; 

(ii) Queue length of the other direction; it is 
irritating when you have to wait in front 
of a red light, but it is even more irritat- 

r ing when there are no (or few) cars pass- 
ing from the other direction. So the queue 
lengths of the other direction and the irri- 
tation per car are inversely proportional; 
(hi) The distance of the car to the traffic light; 
If we assume that the distances between 
two cars in the queue is constant, the dis- 
tance to the traffic light can be written in 
terms of the position k of the car in the 
queue. We have in mind here two concep- 
tually different scenarios: 

— The closer you are to the traffic light, 
the higher the irritation because if the 
light had been green a little longer, you 
were able to pass the roadblock (we call 
this Case I). In this case, the position 
k in the queue and the irritation per car 
are inversely proportional. We use the 
function f(k) = Ci/k, where Cj > is 
a dimensionless constant (hence f{k) is 
dimensionless); 

— The further away you are from the traf- 
fic light at the crucial moment, the higher 
your irritation (there are still so many 
cars in front of you) . We call this Case 
II. In this case, the position k in the 
queue and the irritation per car are di- 
rectly proportional. We use the func- 
tion f(k) — Cuk, where C/j > is a 
dimensionless constant (hence f(k) is 
dimensionless). 

Assuming that (i), (ii) and (hi) affect irrita- 
tion independently, gives us the following ex- 
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pressions for the irritation per car in Direction 

1 (i (1) (fc)) and 2 (i^(k)): 



f(k) 



N 2 {nT- 



;( 2 ) 



(A, TO := 



-Ti)- 
/(*) 



1 



JVi((n+l)T-T) + l 



(T - Ti), 



•2r). 
(5) 

This irritation per car of each direction has units 
s • car -1 because f(k) is dimensionless, iVj has 
dimension car and T, T\ have dimension seconds 
for j = 1,2. Taking the sum of the irritation per 
car (in both directions) defines the total irrita- 
tion at the crucial moment. But for each direc- 
tion, there is only one crucial moment. At all the 
other moments, the light stays red or switches 
from red to green. At these moments, we as- 
sume that there is no extra irritation per car, 
only the irritation that arises from the cars that 
have to wait because they have arrived during a 
red phase. This gives the following expression: 



NrinT+T^ 

= £ ^ 

k=i 



: (1) (fc,Ti) + ^ai(T- 



(6) 

which is the sum of the irritation per car of di- 
rection 1 at time nT-\-T\ plus the irritation that 
occurs from the waiting time of cars that arrive 
in the red phase which lasted T — T\ seconds. 
Note that it is needed to multiply the second 
term with a constant C\ (with units car~ 2 ) to 
make sure that the terms entering (6) (and (7)) 
have the same units. Similarly, we have: 



JV 2 ((n+l)T- 

4 2) (Ti) := X) 

k=l 



,'(2) 



(k,T 1 )+ < fa 2 (T 1 +2T) 



(7) 

which is the sum of the irritation per car over 
all the cars in queue 2 plus the irritation that 
occurs from the waiting time of cars that arrive 
during the red phases of direction 2. Again, C2 
is a constant with units car~ 2 to make sure all 
the terms in (6) have the same units. The total 
irritation in cycle n is then defined as the sum 
of {6| and 0: 



Tift) + 



(8) 



This total irritation has units in s • car -1 . The 
way the irritation is defined makes it very com- 
plicated to minimize this and get a simple ana- 
lytical expression for T\. 

In the next section, we investigate numeri- 
cally the effect of minimizing the total waiting 



time and the total irritation TI using real traffic 
data. 



3 Application: Rush hour in Hillegom 

In this section we use the empirical data that we 
received from Peter Veeke, Traffic and Trans- 
portation Advisory Consultant for Breijn B.V. 
The data contains traffic intensities (in cars per 
hour) during rush hour for some main streets 
in Hillegom, The Netherlands. The traffic in- 
tensities in cars per hour can be converted into 
arrival rates in cars per second. Here, we only 
consider one main street in Hillegom: the Leid- 
sestraat which lies between Singel and Olympia. 
This part of the Leidsestraat is a two-way street. 
The corresponding the arrival rates are given in 
Table 1. We note that during the morning rush 
hour, direction N — > S has a larger arrival rate, 
while in the evening, direction S —> N has a 
larger arrival rate. 



Table 1. Arrival rates for Leidsestraat 



Direction 


Arrival rate 


Arrival rate 




(morning) 


(evening) 


S N 


0.190 


0.264 


N -> S 


0.302 


0.176 



At this point, we assume that there is a road- 
block in the Leidsestraat between Singel and 
Olympia. Given only the arrival rates in Ta- 
ble 1, we cannot determine any values for T\ 
yet. Therefore, we have to make some assump- 
tions about the corresponding passing rates and 
about the cycle and transit period. Because the 
Leidsestraat is a major street and the traffic in- 
tensities are high during rush hour, we assume 
heavy traffic. We take the cycle period equal to 
30 seconds, and the transit period equal to 5 
seconds. We have not received any data about 
the passing rates. Consequently, we need to es- 
timate the passing rates. We assume both di- 
rection have the same passing rates, resulting 
in f3 = j3\ = /?2- Doing the same type of analy- 
sis as in pi, we obtain the following inequalities 
belonging to a heavy traffic scenario: 



(a 2 + 3ai)/?T - 4f3 2 r > 2f3 2 T 



(9) 
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(ai + 3a 2 )f3T - 4/3 2 r > 2/3 2 T (10) 

Hence, our passing rate is typically varying be- 
tween to 0.2. Larger values would simply mean 
that the traffic scenario cannot be considered as 
heavy traffic anymore. Plotting the total wait- 
ing time as a function of the value for T\ and 
the passing rate /3 = fa = fa give us Figure [2] 
and 03 




Fig. 2. 3D plot of the total waiting time against 
the different green times and passing rates for the 
Leidsestraat during morning rush hour. 




Fig. 3. 3D plot of the total waiting time against 
the different green times and passing rates for the 
Leidsestraat during evening rush hour. 



Obviously, we see that the largest value for (3 
results in the smallest total waiting time. But it 
is not obvious that the largest value for j3 also 
results in the smallest irritation, because of the 
construction of the irritation per car. So, we plot 
the irritation in case I as a function of /3 and T\ 
which results in Figure [4] for the morning and [5] 
for the evening rush hour scenario. For the ir- 
ritation in case II we obtain Figure [6] and [7] for 



morning and evening rush hour, respectively. 



00 




Fig. 4. 3D plot of the total irritation in Case I 
against the different green times and passing rates 
for the Leidsestraat during morning rush hour. 




TVs) '' 



Fig. 5. 3D plot of the total irritation in Case I 
against the different green times and passing rates 
for the Leidsestraat during evening rush hour. 

We see in all figures that j3 = 0.2 results in 
a minimal total waiting time. Therefore, we fix 
the passing rates to 0.2. We note that if we fix j3 
and let a increase, we obtain a larger total wait- 
ing time and irritation. This is due to fact that 
increasing a is the opposite effect of increasing 
fa Most likely, cars queues grow unboundently 
and the traffic enters the unstable regime. 
Now, we are able to calculate the values for the 
green-times T\. Because the total waiting time 
and the irritation have different dimensions and 
magnitudes, plotting them against the T\ will 
not gives us a nice figure. If we plot the total 
waiting time and the irritation in Case I and II 
divided by their value at T\ = 0, we get func- 
tions which start at the point (0, 1) that we use 
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Relative value 



150 000 
TI case II 
100 000 



Fig. 6. 3D plot of the total irritation in Case II 
against the different green times and passing rates 
for the Leidsestraat during morning rush hour. 




Fig. 7. 3D plot of the total irritation in Case II 
against the different green times and passing rates 
for the Leidsestraat during evening rush hour. 



to compare the locations of the minima; see Fig- 
ure [8] and [9] Figure [8] shows that the total wait- 
ing time during morning rush hours is minimal 
for Ti = 5.8. Minimizing the irritation in Cases 
I and II gives us T\ = 5.4 and Ti = 3.7 respec- 
tively. For the evening rush hour, minimizing 
the total waiting time and the irritation in Case 
I and II, leads to T x = 13.3, T x = 14.0 and 
Ti = 15.0 respectively, see Figure [9] They are 
different clearing strategies - a priori it is not 
clear cut which one is the best. 



4 Discussion 

We used a discrete model to describe the cars 
traffic in the presence of a traffic lights-controlled 
roadblock under heavy traffic flow conditions. 
Based on our model, we get the expressions ([9| 
and ( 10 ) which are equivalent to those obtained 
in [5j for the two one-way crossroads example. 
This indicates that the two one-way crossroad 
and the roadblock have a similar mathematical 




Fig. 8. Irritation and total waiting time for Leid- 
sestraat during morning rush hour, where 
T/(Ti) . 



1 is given by 



TJ(0) 



in Case I vs. Ti, 



r rj( r r \ 

2 is given by Tj^jz^ m Case II vs. T\ and 

_. . , TWT{T{) 

3 18 glVenby TWT(0) yS - Tl - 




Fig. 9. Irritation and total waiting time for Leid- 
sestraat during evening rush hour, where 
r pj( r p \ 

1 is given by — ^7— ^ in Case I vs. T\. 

r pj( r p \ 

2 is given by jij(q) m Case II vs. Ti and 

TWTiTi) 

3 is given by — — — — — vs. T\. 

6 y TWT(0) 



structure. 

For constant arrival and passing rates we obtain 
Q as explicit expression for the green times at 
which the total waiting time is minimized. As 
main novelty, we introduce two models for the 
driver's irritation which finally lead to two clear- 
ing policies. The driver's irritation depends here 
linearly on the waiting time. At the first sight, 
this is reasonable because the longer the driver 
has to wait, the more irritated he or she will be. 
Such linear vs. nonlinear dependencies need to 
be further explored at a later stage, involving 
the input of group psychologist^] Here we as- 
sumed that the irritation per driver depends on 
the position of that driver in the queue at the 



2 We hope that group psychology is able to bring 
more behavioral insight in the modeling of |5]). 
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moment the light switches from green to red. 
This can be explained in two distinct ways: The 
closer you are to the traffic light, the more ir- 
ritated you become - it was almost your turn 
to pass and now you have to wait again (Case 
I)! Alternatively, the closer you are to the traf- 
fic light, the less irritated you become - now, 
it is almost your turn to pass. The further you 
are from the traffic light the more irritated you 
are because you are wondering why the light 
switches (Case II). These modeling choices are 
based on how we experience irritation at traffic 
lights. There are many ways to model driver's 
irritation, but which is the right one? More re- 
search is needed for a good understanding of 
driver's irritation. Note that the numerical re- 
sults reported in Section [3] indicate that mini- 
mizing the total waiting time looks quite similar 
to minimizing the irritation in Case I; the irrita- 
tion in Case II is leading to somewhat different 
results. 

We applied our model to the traffic data re- 
ceived from Hillegom and found different val- 
ues for the green times in the heavy traffic state 
when minimizing the total waiting time and the 
irritation. Therefore our model raises the natu- 
ral question: 

What is the best strategy: maximize traffic 
throughput or minimize total driver's irritation? 

It is worth mentioning that, unlike |8 (the 
Kumar- Seidman network for production lanes) 
and 10 , e.g., our model is unable to preserve 
bounded waiting queues at roadblocks, and there- 
fore can be classified as unstable (cf. [9], e.g.). 
We hope to be able to repair this within a modi- 
fied framework: reformulate the roadblock prob- 
lem in terms of a discrete particle system, where 
the irritation per particle (car) will be designed 
to depend also on the neighbors positions and 
velocities, which is at the moment not sufficiently 
well captured by ([5|; see 4 , e.g. for a concrete 
possible direction. Such a modeling strategy is 
remotely resembling the stochastic dynamics of 
active Brownian particles. As a next step, we in- 
tend also to apply this modeling framework to 
optimizing for load levels within capacity and 
will continue with further exploration of the de- 
pendence of optimization on modeling of the 
concept of "driver's irritation". 
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